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U: Abstract 

q 

' Motivated by questions in robust control and switched linear dynamical systems, we consider the 

problem checking whether all convex combinations of k matrices in /J"^" are stable. In particular, we 

are interested whether there exist algorithms which can solve this problem in time polynomial in n and 

Ti^lj- ■ k. We show that if fc = [n'^] for any fixed real d > 0, then the problem is NP-hard, meaning that no 

> 



polynomial-time algorithm in n exists provided that P ^ NP, a widely believed conjecture in computer 
science. On the other hand, when fc is a constant independent of n, then it is known that the problem 
may be solved in polynomial time in n. Using these results and the method of measurable switching 



I. Introduction 



\^ i rules, we prove our main statement: verifying the absolute asymptotic stability of a continuous-time 

' switched linear system with more than n'^ matrices Ai e ^"X" satisfying ^ + Aj is NP-hard. 

a 

\ Let {y4i, . . . , Ak] be a finite set of real n x n matrices. We define the corresponding matrix 

c3 ■ polytope A as those matrices which can be written as J2i=i O-iAi for some nonnegative real 
numbers ai, . . . ,0;^ adding up to 1. We are concerned with the following decision problem, 
which we will call the [k, n)-POLYTOPE-STABILITY problem: given k rational nxn matrices 
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{Ai, Afc} to decide whether every matrix in A is Hurwitz stable, i.e. has eigenvalues with 
negative real parts. 

The (/c, n)-POLYTOPE-STABILITY problem, and some of its natural generalizations, have 
been considered before in the control theory literature (see [1], [11], [7], [24], [29], [10], [2], 
[18], [16], [4], [27], [9], [8],[21], [25]). Notable results include the solution of the A; = 2 case 
in [1],[11], [24], also known as the stability testing of affine representations; a Lyapunov-search 
algorithm in [10]; and a recent approach based on LMI relaxations [16], [9], [27]. 

Our first result is that when there are \n'^~\ different nxn rational matrices A^, where d is any 
positive real number (in the above notation, the ([n'^],n)-POLYTOPE-STABILITY problem) 
the problem of deciding whether there exists an unstable matrix in A is NP-hard. 

On the other hand, in many circumstances where the {k, n)-POLYTOPE-STABILITY problem 
appears, k is constant that does not vary with n. Consider, for example, testing the stability of 
an nx n interval matrix where all but r of the entries are known precisely. It is easy to see that 
this is a special case of the poly tope stability problem with k = T' . Moreover, if r is fixed, but 
n is allowed to vary, we end up with a poly tope stability problem with fixed k. 

We note that the {k, n)-POLYTOPE-STABILITY problem can be solved in vP^''^ elementary 
operations by reducing the problem to deciding whether a certain multivariate polynomial does 
not have real roots, which can in turn be solved using quantifier elimination algorithms - see [15] 
for a detailed writeup (the reduction to a root localization problem for a multivariate polynomial 
was first done in [16]). If the number k is fixed, this gives a deterministic polynomial-time 
algorithm. 

The second and main subject of our paper is the {k, n) - Continuous Time Absolute Switching 
Stability Problem ((/c, n) — CTASS problem): given k rational matrices ^i, . . . , in i?"^" to 
decide whether there exists a norm 1 1 • 1 1 in i?" and a > such that the induced operator norms 
satisfy the inequalities: 

\\exp{Ait)\\ < e""* : a > 0, i > 0, l<i<k. (1) 

We consider matrices that satisfy the further restriction of being non-strict contractions in the 
two-norm: 

\\exp{Ait)\\2 <l,t>0,l <i < k, 
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which is equivalent to requiring 



^ A + A*. 



We show that, even in this case, checking the condition of Eq. ([T]) with more than n''- matrices 
Ai G i?"^", where d is any positive real number, is also NP-hard. As far as we know, this is the 
first hardness result for continuous-time switched linear systems. 

We note that it is well known that the stability of the polytope A is necessary, but not in general 
sufficient, for the absolute switching stability condition of Eq. ([U) to hold. Luckily, these two 
conditions are equivalent for the "gadgets" used in our proof of NP-HARDNESS of polytope 
stability. 

We stress that our result says nothing about the solvability of concrete, finite-size instances 
for the problem - for instance, this paper has no new implications for testing the stability of all 
convex combinations of 3 matrices in i?^^^. Rather our result implies that if P 7^ NP then any 
algorithm for solving this problem with n'^ matrices n x n matrices, for any d > 0, will have a 
worst-case operations count which grows faster than any polynomial in n. 

This provides an explanation why many approaches to this problem fail. Despite an extensive 
literature devoted to this problem cited above, no polynomial time algorithms are known. See, 
in particular, [7], for proofs that various intuitive approaches fail. Our result implies that in 
fact any polynomial time algorithm for this problem would immediately disprove the P 7^ NP 
conjecture. 



In this section, we consider the computational complexity of deciding whether every matrix 
in the set A (defined by Eq. (??)) is stable. We will show that this problem is NP-hard through 
a reduction from the maximum clique problem, which is known to be NP-complete [17]. The 
details of this reduction are described below. 

Theorem 1: {n, 2n + 2)-P0LYT0PE-STABILITY is NP-hard. 



Notice that the well-known interval stability proble corresponds to (/c, m)-POLYTOPE- 
STABILITY with k exponential in m. In other words NP-hardness of interval stability, shown 
in [20] and [23], does not imply NP-hardness of (n, 2n + 2)-P0LYT0PE-STABILITY. 

'The interval stability problem is to determine, given numbers {Sij, aij}i,j=i...m, whether every matrix A satisfying Aij £ 
[a^j, Qij] is stable. 



II. NP-Hardness of Stability Testing of Matrix Polytopes 
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We first give a series of definition and lemmas before proving Theorem [T] Given k rational 
matrices Ai, . . . ,Ak in i?"^", we will refer to the problem of deciding whether there exists a 
singular matrix in the associated polytope A as the (A;, n)-POLYTOPE-NONSINGULARITY 
problem. 

Lemma 2: There is a polynomial-time reduction from the {k, n)-POLYTOPE-NONSINGULARITY 
problem to the {k, 2n)-P0LYT0PE-STABILITY problem 
Proof: Given a square matrix A E R,"X" define B as 



B 



n 4^ 



A Inxn 

We claim B is Hurwitz if and only if A is nonsingular. This follows because spectrum of B can 



be written as cr(_B) = Ui<i<„{|(— 1 ± y 1 — where Si{A) < ... < s„(A) are singular 

values of the matrix A. 

Suppose we are given k nx n matrices A^, and we want to decide whether the set A defined 
by Eq. (??) contains a nonsingular matrix. Define 



Since 



V 



n A' 

'-'nxn ^ 



T 

i 

Ai Inx 



^oiiBi 



V 



Onxra 

Si (^iAi Inxn 

when Eitti = 1, it follows by the previous item that testing POLYTOPE-NONSINGULARITY 
with the set A is the same as testing POLYTOPE-STABILITY on the set 

B = {B \ B = OiBi, ^ = 1, ftj > for all i}. 

i i 

However, note that the construction has doubled the dimension, since the matrices Bi belong to 
^2nx2n jj^j^ concludcs the proof that (A;, n)-POLYTOPE-NONSINGULARITY can be reduced 
to {k, 2n) -POLYTOPE-STABILITY. ■ 
Consider the problem of deciding whether there exists a nonnegative vector p in i?" whose 
components sum to 1 such that p^Mp = 1 for an arbitrary invertible matrix M. We will consider 
to be the input to this problem. We will refer to the problem as the n-QUADRATIC- 
THRESHOLD problem. 
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Lemma 3: There is a polynomial-time reduction from the n-QUADRATIC-THRESHOLD 
problem to the {n, n + l)-POLYTOPE-NONSINGULARITY problem. 
Proof: 
Define 

;,(n) ^ { 

' \ ef 1 

where is the column vector with 1 in the i'th entry and zeros elsewhere. Define <S„ = {p e 

R"" I EiPi = 1, Pi > Vi} and let 

i 

In other words, X is the set of matrices of the form 

[ P 1 J' 

with p e »S„. By the Schur complement formula such a matrix is singular if and only if p^Mp ~ 
1. Thus given an invertible matrix M, we can solve the n-QUADRATIC-THRESHOLD problem 
by solving an instance of the n+l-POLYTOPE-NONSINGULARITY problem with the polytope 
X. U 

Let G = {V,E) be an undirected graph. A subset of the vertices C is called a clique if 
Vi,V2 E C implies {vi,V2) E E. The n-MAX-CLIQUE problem is the problem of determining 
the size of the largest clique, denoted by uj{G), in an undirected graph G on n vertices. 

Lemma 4: There is a polynomial-time reduction from the n-MAX-CLIQUE problem to the 
n-QUADRATIC-THRESHOLD problem 
Proof: 

1. It is known that [19]: 

1 J. 

1 -— = maxp Mp. (2) 

uj{G) peSn 

where M is the adjacency matrix of the graph G. 

2. Because the QUADRATIC-THRESHOLD problem is defined only for nonsingular matrices, 

for our reduction to work we will need to modify M to insure its nonsingularity. To this end, 
we consider the matrices Mi = M + ^^s^/ for i = 1, . . . , n + 1. At least one Mj must be 
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nonsingular, because M cannot have n + 1 eigenvalues. We find a nonsingular Mi (tliis can be 
done in polynomial time with Gaussian elimination for each i = l,...,n + l). Let us denote 
this nonsingular Mi by Mi*. 

In the proof below, we will threshold the form p'^Mi*p; recall from Lemma [3] that this requires 
the computation of Af^T^. This involves a polynomial number of computations in n, and the bit- 
sizes remain polynomial as well; for a proof see Corollary 3.2a of [26]. 
3. We have that for p E Sn, 

p^M,,p = p^Mp + P^i 
Because pi G [0, 1] and I^iPi = 1 imply that Y^iPi < 1^ we have: 

p^Mp < p^Mi^p < p^Mp + — (3) 

+ I* 



It follows that 



1 ^ 11 

< maxp^ Mi,p < 1 — + — (4) 



4. Because the optimal solution of Eq. (O is 1 — l/uj(G), and uj(G) is an integer between 1 and 
n, this optimal solution must be in the set S = {0, 1/2, 2/3, 3/4, . . . , 1 — l/n}. Because the gap 
between the elements of S is less than l/n^, and consequently less than l/(n^ + i*), it follows 
from Eq. ^ that the largest element of S smaller than maxpg^^ p^Mj.p must be 1 — ^j^, the 
solution of the MAX-CLIQUE problem. Let this element be called k*; then, it follows that 

^P^Mi*p > 1. (5) 

for some p E Sn, and k* is the largest element of S with this property. 

For each k E S, the existence of a p E Sn satisfying Eq. ^ can, due to the invertibility 
of M*, be decided by evaluating p^Mi*p at an arbitrary p E Sn, followed up with a call 
to the QUADRATIC-THRESHOLD problem. This is the reduction from MAX-CLIQUE to 
QUADRATIC-THRESHOLD. ■ 

Proof: [Proof of Theorem [T]|: Lemmas [2l[3l|4] provide a reduction from the MAX-CLIQUE 
problem to the POLYTOPE-STABILITY problem. The size of the problem goes from n in 
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the QUADRATIC-THRESHOLD problem to (n,n + 1) after the reduction to POLYTOPE- 
NONSINGULARITY; and from {n, n + 1) to {n, 2(n + 1)) in the reduction from POLYTOPE- 
NONSINGULARITY to POLYTOPE-STABILITY. Since MAX-CLIQUE is known to be NP- 
complete [17], it follows that (n, 2n + 2)-POLYTOPE-STABILITY is NP-hard. ■ 
Remark: Note that the matrices Ai created in our reduction have entries whose bit-size are 
polynomial in n. 
Remark: 

1) We note that our results easily imply the NP-HARDNESS of POLYTOPE-STABILITY 
with In^l extreme points, for any real d > (here [x] refers to the smallest integer 
which is at least x). Indeed, the [n'^] -MAX-CLIQUE problem remains NP-COMPLETE, 
and therefore the ( [n'*] , 2 [n'*] +2) -POLYTOPE-STABILITY remains NP-Hard. However, 
we clearly do not make the problem any easier by increasing the dimension, so that the 
( \n'^] , 2n + 2) problem is NP-Hard. 

2) We have remarked that the (A;, n)-POLYTOPE STABILITY problem may be solved in 
polynomial time if k is upper bounded by a constant. We may ask about the reverse 
question: what happens when n is upper bounded by a constant? Is the problem solvable 
in polynomial time as a function of A;? 

The answer is yes. Caratheodory's theorem implies that any matrix in A may be expressed 
as a convex combination of + 1 matrices. Thus we reduce the problem of checking {k, n) 
polytope stability to the problem of checking {^^2_^^ different {ri^ + 1, t^) polytope stability 
problems. When n is upper bounded, checking {n^ n) -POLYTOPE- STABILITY (say 
by computing the determinant explicitly and using quantifier elimination) takes a constant 
number of operations, so the number of operations grows as („2^i), which, when n is 
upper bounded, is polynomial in k. 

3) Essentialy, our construction boils down to the next determinantal representation: 

<5(Pi,-,Pn) =det(Ao+ ^ PiA) (6) 

l<i<n 

Such representations exists for any polynomial Q{pi, ...,Pn) [28]. In our case 

g(pi,...,p„) =< Mp,p> -L 
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III. NP-HARDNESS OF Checking Continuous-Time Absolute Switching 



Stability 



In this section, we will show that checking the absolute switching stability of a class of 
continuous-time linear switched systems is NP-hard. Recall that the (k, n) continuous time 
switching stability problem is: given k rational matrices Ai, . . . , in i?"^", to decide whether 
there exists a norm || ■ || in i?" and a > such that the induced operator norms satisfy the 
inequalities: 



We consider a subcase of the problem where the matrices Ai satisfy the (nonstrict) Lyapunov 
inequalities ^ + , \ < i < k. In assuming this, we only make the problem easier, 
as we assume that the matrices have a nice geometrical structure; indeed, the previous 
requirement corresponds to requiring that solutions of the equation x{t) = Aix{t),i = 1, . . . ,k 
with measurable switching rules, are nonincreasing in the 2-norm. 

Nevertheless, we will show that testing continuous time stability is NP-hard already in this 
case. The following lemma - together with Theorem [T]- proves this. As far as we know it is the 
first hardness result in the area of continuous time absolute switching stability . 

Lemma 5: Consider the following 2n x 2n matrices 1 < i < A; < oo: 



Then there exists a norm 1 1 ■ 1 1 in i?^" and a > such that the induced operator norms satisfy 
the following inequalities: 



if and only if all matrices in the convex hull A are nonsingular. 

Before proving the lemma, we need the following auxiliary claim. Consider the following 



exp{Ait)\\ < e""* ■.t>0,l<i<k. 



(7) 




\\exp{Bit)\\ < e""* ■.t>0,l<i<k, 



(8) 
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family of "differential" equationj^: 

k 
i=l 

with initial condition satisfying ||a:;(0)||2 = 1. Here ^ pi{t) . . . Pk{t) ^ a Lebesgue-measurable 
vector function whose range is a subset of Sk- Since the Y^iPiBi lies in a bounded set of 
matrices, the above equation has a unique Lipschitz solution. Since ^ 5^ + Bj , it follows that 

||x(t)||2 < 1 for t > 0. 

Claim: Assume that there is no induced norm satisfying Eq. Then, there exists a measurable 
vector function ^ Pi{t) . . . Pk{t) ^ whose range is a subset of Sk, and a vector x(0) with 
||x(0)||2 = 1 such that the solution of 

iit) = [Ep.(t)i?.]x(t), (9) 

i=l 

with initial condition a;(0) satisfies ||2;(1)||2 = 1. 

Proof of claim: We will prove the contrapositive of the claim. It can be seen that the set 
of all possible values of x{t) at time t = 1 produced by choices of p(t),x(0) which satisfy 
our assumptions is compact; see Theorem 4.7 in [13] for details of the proof. So, suppose the 
conclusion is not true, by compactness this means there exists 5 > such that for every p, x{0) 
satisfying our assumptions, we have that ||x(l)||2 < I — 5. Thus there exists an e > such that 
the system 

k 

= [J2P^it)iB^ + eI)]x{t), (10) 
i=l 

has the same property( i.e. there exists 6 with ||x(l)||2 < 1 — 5 for all suitable choices of 
x(0),p(t)). 

We will define a norm such that ||e^'*|| < e~"* for some a > and all t > 0, thus violating 
Eq. ([8]). Define a norm on i?" as follows. For q E -R", 

||g||„ = sup sup l|a;(ti)||2 

'^^^ measurable p{z),ze[o,u] with range in 

^Strictly speaking, they ought to be viewed as integral equations 



x(t)=x{0)+ / {["^p^{r)B,]x{r))dr. 
-^0 ,=1 
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where x{-) is a solution to Eq. (flOl) with a;(0) = q. 
This norm induces a norm on i?"^": 

\\Q\\nxn= sup IIQxIln 
3::6i?"||a::||„ = l 

However, 

< IItII 

since premultiplication by e^'^"^ corresponds to simply requiring that p{z) = Cj for the first t 
time units. Therefore, 

\\pBa\\ < p-^i- 

I 1"^ I Inxn _ ^ 

which proves the claim. 

Proof: [Proof of Lemma [Sj] We remark that the argument is very similar to Theorem 4.7 
in [13] and Corollary 2.8 in [14]. 

First, we show the "only if" part, that is, assuming Eq. ([8]), all the matrices in the convex 
hull A must be nonsingular. Indeed, suppose not; suppose there exists a vector p such that 
Yli=iPiBi is not stable, with J2i=iPi = 1 and all Pi > 0. Let A be an eigenvalue of J2iPiBi with 
nonnegative real part; then, will be an eigenvalue of e^*^'^'. Thus, e^»^'^' has an eigenvalue 
of magnitude at least 1, so that 

||eE.P»^«|| > 1, (11) 

where || ■ || is the same norm as in Eq. ([8]). However, by the Baker-Campbell-Hausdorff formula, 
for all e > 0, there exists m large enough so that 

||(eE.=iP«s»)|| < (l + e)||(J]e(^/'")P'^0™|| 

i=l 

Since | |e(i/™)p»^»| |™ < e^'^^^ by Eq. ([8]), we have that large enough m 

||eEtiP»s«|| < (i + e)e-'^E.P' = (1 + e)e-'' < 1, (12) 



where we pick e small enough for the last inequality to hold. Equations (fTTI) and (fT2l) are in 
contradiction. We conclude that all the matrices in B are indeed Hurwitz. By Lemma [2l this 
implies all matrices in A are nonsingular. This proves the "only if" part. 
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Next, we show the "if" part. Let x(0),p(t) be such that the conclusion of the above claim is 
satisfied. The corresponding curve x{t) : < t < 1, satisfying Q is Lipschitz; ||x(t)||2 = 1 for 
all < t < 1. It follows from the structure of matrices Bi that if C E B then for all non-zero 
vectors x the inner product {Cx,x) < and {Cx,x) = if and only if a; G i?" © 0„. Here 
-R" © On is the subspace spanned by {ei, . . . , e„}. If for some r G [0, 1) the vector x(r) does 
not belong to i?" © 0„ then it also holds by continuity of the curve x{t) in a sufficiently small 
neighbourhood [r, r + e]. This implies that {x{t),Cx(t + 5)) < —b < for some b > and 
all < (5 < e and C G i3. Define the vector v{e) = J^^^{[Ei=iPiis)Bi]x{s))ds. It is clear that 
||?^(e)||2 < Ke for some constant K and that (x(r),i;(e)) < —be. Therefore 

||x(r + e)||^ = ||x(r)||^ + 2 < x(T),t;(e) > + 

\He)\\2 
< l-2be + K'^e^. 

We get for small enough e the inequality ||2;(r + e)||2 < 1. We conclude that x{t) G ® On '■ 
0<t<l. 

This gives that jQ{J2i<i<kPi{'^)^i)^i{''')dT — for all < t < 1, where Xi(t) is the vector 
formed by the first n components of x(r). As the Lebesgue measurable vector function on 
[0,1] ^ p^(t) ... pk(t) ^ e Sk is bounded thus it follows that (J2i<i<kPi{^)A)^i{'^) = 
up to measure zero. Since the last n components of x(r) are zero, and ||x(t)||2 = 1, we have 
\\xi{t)\\2 = 1,0 < t < 1. Therefore there must exist a singular matrix in A. ■ 

Remark: In the discrete time case, it is known that given two nx n rational matrices A, B it is 
NP-hard to check if there exists a norm ||.|| in i?" such that the induced norms \ \A\\, \ \B\ \ < 1 
[5]. On the other hand, it had been observed in [12], that a slight modification of a construction 
in [5] gives a direct proof of the following statement: given two n x n rational matrices A, B 
with nonnegative entries and ||v4||ij, \ \B\\i-^ < 1 it is NP-HARD to check if there exists a norm 
1 1 . 1 1 in i?" such that the induced norms \ \A\\, \ \B\ \ < 1. However, the continuous-time counterpart 
of this last problem is "easy"(see Theorem 2.1 in [14]). Based on this, it is unclear whether it is 
possible to modify the constructions from [5], [12] to handle the continious time case. 
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